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Abstract 


bJT) Strongly-coupled quantum dot molecules (QDMs) are widely deployed in the design of a variety of optoelectronic, 
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photovoltaic, and quantum information devices. An efficient and optimized performance of these devices demands 
engineering of the electronic and optical properties of the underlying QDMs. The application of electric fields 


(N 


offers a knob to realise such control over the QDM characteristics for a desired device operation. We perform 
multi-million-atom atomistic tight-binding calculations to study the influence of electric fields on the electron and 
hole wave function confinements and symmetries, the ground-state transition energies, the band-gap wavelengths, 
g and the optical transition modes. The electrical fields both parallel {Ep) and anti-parallel (Ea) to the growth 
direction are investigated to provide a comprehensive guide on the understanding of the electric field effects. The 
S strain-induced asymmetry of the hybridized electron states is found to be weak and can be balanced by applying 
^ a small Ea electric field, of the order of 1 KV/cm. The strong interdot couplings completely break down at large 

o 

O electric fields, leading to single QD states confined at the opposite edges of the QDM. This mimics a transformation 
I from a type-I band structure to a type-II band structure for the QDMs, which is a critical requirement for the 

^ design of intermediate-band solar cells (IBSC). The analysis of the field-dependent ground-state transition energies 

00 

reveal that the QDM can be operated both as a high dipole moment device by applying large electric fields and as 

in 

a high polarizibility device under the application of small electric field magnitudes. The quantum confined Stark 

00 

effect (QCSE) red shifts the band-gap wavelength to 1.3 ym at the 15 KV/cm electric field; however the reduced 

in 
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electron-hole wave function overlaps lead to a decrease in the interband optical transition strengths by roughly 
three orders of magnitude. The study of the polarisation-resolved optical modes indicate the benefit of applying 


^ small electric fields, which lead to an isotropic polarisation response, a desirable property for the semiconductor 


a 


optical amplifiers (SOAs). 
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1 — Introduction 


Just as the quantum dots (QDs) are widely known as “artificial atoms” due to their discrete energy 
spectra and the three-dimensional confinement of the charge carriers, the strongly-coupled quantum dots 
are often referred as “artificial molecules” due to their hybridized bonding and anti-bonding type electronic 
states, closely resembling to that of the molecules made up of atoms. Self-assembled quantum dot molecules 
(QDMs) made up of III-V materials are one of the most popular types of nanostructures used in a variety 
of devices for optoelectronic [1-7], photovoltaic [8-12], and quantum information technologies [13, 14], The 
QDMs typically grow in the form of vertical stacks due to the presence of strain, which stems from the 
lattice mismatch of the substrate and the QD material. One such example is of InAs/GaAs QDMs where 
about 7% lattice mismatch between the InAs and GaAs materials results in the formation of nearly perfectly 
aligned, vertically stacked QD layers of InAs when epitaxially grown on top of a GaAs substrate [15, 16]. 
Similar stacking of the QD layers forming QDMs has been demonstrated in other material systems such 
as InAs/InGaAsP/InP [17, 18], InGaAs/GaAs [19, 20], GaSb/GaAs [21], and InGaN/GaN [22, 23], etc. 
The stacking of multiple QD layers offers many useful benefits from the perspective of devices, such as 
increased number of photons emitted or absorbed per unit area at the QD energy levels for photo-diodes 
and lasers, a wider range of light absorption wavelengths for intermediate-band solar cells (IBSGs), access 
to long wavelengths for mid and far-infrared photo-detectors, and polarisation insensitivity of the gain for 
semiconductor optical amplifiers (SOA’s), etc. The performance of these devices is critically related to the 
confinement of charge carriers inside the QDM. The spatial distributions and symmetries of the confined 
electron and hole wave functions control interband/intraband optical emissions and absorptions. Therefore 
it is of crucial significance to properly understand the electronic structure of the QDMs in order to fully 
exploit their promising properties for the design and optimization of devices. 

The electronic structure of self-assembled QDMs is strongly influenced by their geometry and the strain 
which stems from the lattice mismatch of the constituent materials. Therefore, the past theoretical studies 
of the QDMs have primarily focused on these two parameters to investigate the symmetries and the spatial 
confinements of the charge carriers. It has been experimentally shown that the strain not only plays a 
major role in governing inter-dot couplings leading to the formation of molecular states, but also strongly 
influences the polarisation response of the output photoluminescence (PL) spectra, leading to polarisation 
insensitivity of the output light for the properly designed QDMs [3, 24]. For many applications such 
as photodetectors, photovoltaics (such as IBSGs), and quantum information devices, the QDMs are also 
subject to electric fields, which can significantly influence their electronic structure. Such fields are either 
present internally such as investigated recently for the solar cells [25, 26], or can be applied externally to 
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purposely tune the properties of the devices for a certain operation [14], Although the impact of electric 
helds has been extensively studied in the literature for the single QDs [27-29], the bi-layer QDMs [7, 30-34], 
and for the small QD stacks consisting of three to four QD layers [35], to-date no theoretical understanding 
is available on the electronic and optical properties of the multi-layer QDMs under the influence of electrical 
helds. 

One rather intuitive effect of the electric helds is to Stark shift the energy levels of the QDMs and 
therefore modulate the emission/absorption peak wavelengths. However there are several open questions 
which require rigorous theoretical calculations, such as: (1) How the presence/application of electric helds 
modify the spatial conhnements and symmetries of the electron and hole wave functions?, (2) What is the 
magnitude of the electric held required to break the coupling of quantum dot layers in a QDM?, (3) Does 
the impact of an applied electric held directional on the electronic and optical properties of the QDMs?, 
(4) What is the built-in dipole moment and polarizibility for the QDMs as a function of the electric held 
magnitudes and directions, and (5) What are interband optical transition strengths at large electric helds, 
where the band gap wavelength is red shifted to 1.3 pm. A detailed and reliable theoretical understanding 
of these unknown parameters is critically signihcant for the experimentalists working on the design and 
implementation of cutting-edge technologies based on QDMs. 

This work reports a hrst comprehensive theoretical study of the electronic and optical properties of the 
QDMs under the inhuence of electric helds. The physical insights provided are based on the state-of-the-art 
multi-million atom electronic structure and optical mode simulations, and deliver a reliable understanding 
of the complex physics of the QDMs, which involve subtle interplay between the applied electric helds 
and the strong interdot couplings. In order to answer the aforementioned outstanding questions, we apply 
electric helds along both parallel and anti-parallel directions to the growth direction of a QDM recently 
discussed in the experimental studies [8, 9] and investigate its ehects on the spatial conhnements and 
symmetries of the electron and hole wave functions, band-gap wavelengths, ground state transition energies, 
and polarisation-resolved interband optical modes. Our results show that the application of a small electric 
held anti-parallel to the growth direction is sufficient balance the asymmetric ehect of the underlying strain 
on the electron states. At large electric held magnitudes, the interdot couplings completely break down 
and lead to the formation of single QD electron and holes states, which are conhned at the opposite ends 
of the QDM. This mimics a transformation of the QDM electronic structure from a type-I system to an 
ehectively type-H system, in which the single QD electrons and holes are separated by a large distance, of 
the order of the size of the QDM along the growth direction - a desirable property for the design of the 
IBSCs. The band-gap wavelength is calculated to increase roughly linearly as a function of the applied 
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electric fields, irrespective of its direction for the large electric field magnitudes. However at the low electric 
helds (< 2 KV/cm), the change in the band-gap wavelength is highly dependent on the direction of the 
electric held. An interesting property of the QDMs is revealed by analysing ground state transition energies 
and demonstrating that QDMs exhibit different characteristics under small and large helds: at the small 
helds, the QDMs have large polarizibilities but negligible dipole moments and this behaviour reverses at 
the large helds. The increased spatial separation between the electron and hole wave functions is directly 
rehected in the calculated interband optical transition strengths, which decrease as a function of the applied 
electric helds. At an electric held magnitude of 15 KV/cm required for the band-gap wavelength of 1.3 
/im, the interband optical transition strength is calculated to decrease by nearly three orders of magnitude. 
Therefore very small oscillator gain is expected from the QDMs at the 1.3 pm operating wavelengths. 
However the application of small electric helds, ranging from 2 KV/cm to 5 KV/cm, are predicted to 
be useful for the optical devices demanding isotropic polarisation response such as semiconductor optical 
amplihers. 

The rest of the paper is divided into the following sections: the section 2 briehy layouts the methods used 
in our simulations. The section 3 provides results and related discussions. We start with the section 3.1 by 
briehy reporting the electronic states of the QDM without the application of electric helds. This highlights 
the impact of strain on the electronic structure of the QDM and forms a basis for the remainder of the paper. 
The inhuence of the applied electric helds on the electron and hole wave functions is thoroughly investigated 
in the sections 3.2, 3.3, and 3.4, where we provide details about the held-dependent conhnements and 
symmetries of the electronic states. The section 3.5 predicts the held-induced shifts in the band-gap 
wavelengths and the ground state transition energies. The section 3.6 provides calculations of the interband 
optical transitions for the TE and TM mode polarisations under the application of electric helds. Finally 
we conclude in the section 4. 


2 — Theoretical Framework 

Our investigation of the electric held ehects on the QDM properties is based on multi-million-atom 
simulations of the strain, electronic structure and interband optical transitions. The strain is computed 
from the valence force held (VFF) relaxation of the atoms [36, 37]. A very large simulation domain 
consisting of roughly 43 million atoms is resolved to properly accommodate the ehects of the long-range 
strain helds and boundary ehects. The electron/hole energies and states are obtained by solving a twenty- 
band sp^d^s* tight-binding Hamiltonian [38]. The GaAs box surrounding the QDM have realistic boundary 
conditions as previously reported in Ref. 39, and the surface atoms are passivated according to the published 
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recipe [40]. 

The polarisation-resolved interband optical transition strengths are calculated as follows: first we cal¬ 
culate interband optical transition strengths between the individual electron and hole states by using the 
Fermi’s golden rule, where the squared absolute value of the momentum matrix elements is summed over 
the Kramer spin degenerate states: Te.^hj = |(ei|['n^, H]|hj)p where H is the single particle tight binding 
Hamiltonian in the sp^d^s* basis, e* and hj are the ith and jth electron and hole states respectively, and 

is a vector along the polarisation direction. The polarisation dependent optical modes TEjno], TE[_iio], 
TE[ioo], TE[oio], and TM[ooi] between the lowest conduction band state ei and the highest four valence band 
states hi, h 2 , ha, and h 4 are calculated by rotating the polarisation vector it = {it + lt) cos0sin 9 + ^ cos 6 
to an appropriate direction in the polar coordinates: for the TEpioj-mode: 6 = 90° and 0 = 45°, for the 
TE[_iio]-mode: 6 = 90° and 0 = 135°, for the TE[ioo]-niode: 6 = 90° and 0 = 0°, for the TE[oio]-niode: 
6 = 90° and 0 = 90° and for the TM[ooi]-niode: 6 = 0°. Next each optical transition strength is artificially 
broadened by multiplication with a Gaussian distribution centred at the wavelength of the transition [41]. 
Finally we add all of the four Gaussian functions to calculate the total optical intensity function, f(\). The 
complete expression for the optical intensity function is given by equations 1 and 2 as follows: 


where. 
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To investigate the effect of electric fields, the potential due to a static uniform electric field is added 
in the diagonal of the tight-binding Hamiltonian. We consider electric helds both parallel to the growth 
direction [001] (labelled as Ep) and anti-parallel to the growth direction (labelled as Ea). The magnitude of 
the electric helds is varied between 0 and 45 KV/cm, which is in accordance with the range of the electric 
helds recently investigated in an experimental study for the design of solar cell devices [25]. The atomistic 
simulations are performed using nanoelectronic modeling tool NEMO-3D [42-44]. 

Fig. 1 (a) illustrates the schematic diagram of a QDM in the form of vertically-stacked closely-spaced 
QD layers. The QDM consists of nine QD layers (QDLs) made up of InAs material and is embedded into 
a large GaAs box. The geometry parameters of the QD layers (Shape = Dome, Base diameter = 20 nm. 
Height = 4 nm. Wetting-layer thickness = 0.5 nm, and Spacer thickness = 4.5 nm) are obtained from the 
recent experiments [3, 45], so that our theoretical results remain relevant in the context of the recently 
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FIG. 1. (a) The schematic diagram of the simulated system is shown, illustrating a quantum dot molecule consisting 
of closely-spaced vertically-stacked nine quantum dot layers. Each QD layer is made up of a dome shaped InAs 
QD with the base diameter 20 nm and the height of 4 nm, residing on the top of a 0.5 nm thick InAs wetting 
layer. A large GaAs box surrounds the QDM, where the thicknesses of both the substrate and the cap layer are 
30 nm. The lateral size of the GaAs box is 50 nm. The total number of atoms in the simulation domain is roughly 
43 million. The direction of applied field is also labelled, with the field parallel to the growth axis ([001]) labelled 
as Ep and the field anti-parallel (opposite) to the growth axis is labelled as Ea- (b) The side views of the lowest 
three electron states (el, e2, e3) and the highest four hole states (hi, h2, h3, h4) are shown, indicating the spatial 
confinements of the wave functions. The dashed lines mark the locations of the InAs wetting layers, (c) The top 
views of the lowest three electron states (el, e2, e3) and the highest four hole states (hi, h2, h3, h4) are shown, 
highlighting the symmetries of the wave functions. The dashed circles indicate the boundaries of the QD layers. 


published experiments. The experimental studies [16, 25] have not reported any In-Ga intermixing for 
the QDMs. Therefore we have ignored this effect in our simulations. However, our recent studies [46, 47] 
have indicated that the In-Ga intermixing have minor impact on the conhnements and symmetries of the 
electron and hole wave functions for the single and trilayer QD systems. Gonsequently, we expect that the 
presented physical insights will remain valid in the presence of weak In-Ga intermixings. 


3 — Results and Discussion 

3.1 — A brief overview of the strain effects without the application of electric fields 

Before investigating the impact of electric helds, we briefly describe the conhnement of the electron and 
hole states in the strongly-coupled QDM (as shown by the schematic diagram of Fig. 1 (a)). The physical 
insights presented in this section will serve as a basis to understand the effects of electric helds in the later 
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sections of this paper. The lowest three electron states (el, e2, e3) and the highest four hole states (hi, h2, 
h3, h4) are plotted in Fig. 1, exhibiting two different angles: (b) the side views in the [010]-plane to show 
the conhnement of the states and (c) the top views in the [001]-plane to highlight the symmetries of the 
states. Here a major role is played by the underlying strain arising from the lattice mismatch of the InAs 
and GaAs materials, which strongly impacts both the electron and hole states. The impact of strain on the 
electron states is described in two ways: (1) It strongly couples the QD layers transforming the electron 
states of individual dots into hybridized molecular states spread over multiple QD layers — for example el is 
a bonding-like state and e2 is an anti-bonding-like state and both states are spread over several QDLs. The 
top views of the el, e2, and e3 states further illustrate that all of the three lowest electron states exhibit 
symmetric s-orbital shape. This is in stark contrast to the atomic-like states of the isolated single QD 
layers, where only the lowest el state is of s-orbital symmetry, and the hrst two excited (e2 and e3) states 
exhibit p-orbital symmetries [39]. (2) The second influence of the strain is evident from the asymmetric 
distribution of the electron states over the span of the QDM. For instance, the lowest electron state (el) is 
quite prominently pushed towards the bottom of the QDM. Consequently, the anti-bonding-like e2 electron 
state is pushed towards the top of the QDM. This effect is due to the dome shapes of the QD layers, which 
lead to an asymmetric distribution of strain along the growth direction ([001]) and favours the bottom of 
the QDM for the lowest electron state. 

Contrary to the electron states, which are hybridized over several QD layers, the hole states are much 
more conhned, even in this case of this very strongly-coupled QDM where the inter-layer separation is 
only 1.5 nm. This behaviour of the hole states is a direct consequence of their heavier effective mass, 
which reduces the impact of strain and electronic couplings. The side views of the hole states clearly 
indicate a restriction for the conhnement of the hole states to either a single QD layer or to a couple of 
adjacent QD layers. Looking at the top views in Fig. 1(c), another interesting property of the QDMs is 
observed: the hole states are conhned close to the interfaces and are aligned along the [-110] direction. 
The conhnement of the hole states close to the QD layer boundaries are related to the strong inter-dot 
coupling of the QD layers which implies that the QDM can be considered as a single nanostructure with 
a very large aspect ratio (height/base). A quantum dot with very large aspect ratio has been shown to 
exhibit the formation of heavy-hole pockets at the interfaces which trap the hole states [39, 48], therefore 
similar ehect is observed here for the QDM. The conhnement of the hole-states along the [-110] direction 
is related to the effect of strain, which enhances the symmetry breaking effect and makes the [110] and the 
[-110] directions inequivalence for the conhnement of the electronic states. It has been previously shown 
that the strain favours the [-110] direction as the preferred alignment direction for the electronic states. 


7 
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FIG. 2. Side views in the [010]-plane are shown for the lowest three electron states (el, e2, e3) as a function of the 
electric field (Ep), applied parallel to the growth direction. The magnitudes of the Ep are provided in the units of 
KV/cm. The dashed lines mark the locations of the InAs wetting layers. The lowest electron state el is rapidly 
pushed towards the lowest QDL. The anti-bonding state e2 maintains its character for very large values of the 
applied electric fields (upto 25 KV/cm). At 35 KV/cm, the coupling of the QDLs completely breaks down and 
the lowest three electron states roughly resemble the states of an isolated InAs QD, confined in the QDLi. 


leading to the lower p-state being aligned along the [-110] direction [5, 31]. Therefore all of the hole states 
in the QDM are aligned along this direction. 

From the aforementioned discussion, it is clear that the study and understanding of the electric held 
control of the electronic states in a strongly-coupled multi-layer QDM is much more involved when compared 
to much simpler case of single QDs. This is because in single QDs, electric held only modihes the symmetries 
of the electronic states and its impact on the spatial conhnements is very weak - only a small increase 
between electron-hole separation has been reported earlier in the literature [28]. In contrast, due to the 
large volume of a QDM, the application of the electric helds can lead to much larger separations of the 
electron-hole states. This characteristic of the QDMs may be useful in the design of photovoltaic devices, 
where the reduced electron-hole recombination rates are expected to enhance charge carrier collection at 
the end terminals. Furthermore, the changes in the symmetries of the electronic states would also be quite 
different for a QDM when compared to a single QD, as all of the three electron states are s-type states 
and all of the hole states are interface-conhned states. It should also be noted that the interplay between 
the strain and electric held is expected to be drastically different for the electron and hole states, as they 
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FIG. 3. Side views in the [010]-plane are shown for the lowest three electron states (el, e 2 , e3) as a function of 
the electric field (Ea), applied anti-parallel to the growth direction. The magnitudes of the Ea are provided in 
the units of KV/cm. The dashed lines mark the locations of the InAs wetting layers. The impact of electric field 
pushes the electrons towards the top of the QDM, opposing the impact of strain, which pushes the electron wave 
functions towards the bottom of the QDM. At 1 KV/cm, the effect of the field and the strain is roughly balanced 
and the spatial distributions of the states become symmetric with respect to the center of the QDM (QDL 5 ). At 45 
KV/cm, the interdot coupling of the QDM completely breaks down and the lowest three states roughly resemble 
the states of an isolated InAs QD, confined in the QDLg. 

have very different strengths of conhnement in the QDM. In the next three sections, we will systematically 
investigate the impact of electric fields on the confined electron and hole states, separately considering 
their spatial conhnements and symmetries. 


3.2 — Electric field effect on the electron states 

In this section, we apply nniform electric helds to the QDM and calculate their impact on the lowest 
three electron states. As presented in the last section, the underlying strain held pushes the lowest electron 
state el towards the bottom of the QDM and the effect is in the opposite direction for the second electron 
state e2. Due to this inherent asymmetry of the electron states along the [001] direction, the impact of the 
applied electric held is also expected to be directional: a held parallel to the growth direction (Ep) would 
support the ehect of strain and a held anti-parallel to the growth direction {Ea) will tend to oppose the 
ehect of strain, cancelling its impact on the electronic structure of the QDM. 

Fig. 2 (a-j) shows the lowest three electron states as a function of the electric held magnitude {\Ep\) 
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applied along the growth direction [001], where the magnitude of the held is varied from 0 to 35 KV/cm. 
The electric held pushes the ground state el towards the bottom of the QDM, and therefore supports the 
ehect of strain, which also pushes el towards the bottom of the QDM. The ehect of electric held anti¬ 
parallel to the growth direction Ea is illustrated in Fig. 3 (a-k). In this case, the held pushes the ground 
state el towards the top of the QDM and therefore acts against the ehect of the strain. As evident from 
Fig. 3 (b), the application of a small anti-parallel electric held, \Ea\ = 1 KV/cm, nearly entirely cancels 
the ehect of strain, leading to a symmetric distribution of the lowest electron state along the growth axis 
of the QDM. This property is further highlighted in Fig. 4, where we plot one dimensional (line) cuts of 
the el charge density (IV’eiP) as a function of the distance along the [OOlj-axis through the center of the 
QDM. In the absence of both strain and electric held (shown as blue line), the charge density distribution 
is symmetric with its peak value at the center of the hfth (middle) QD layer of the QDM (QDL 5 ). The 
presence of the strain (in the absence of any electric held) slightly pushes the charge density towards the 
bottom of the QDM (as shown by the red line), with the peak of the charge density now residing in the 
fourth QD layer (QDL 4 ). The application of \Ea\ = 1 KV/cm cancels the ehect of strain, restoring the 
symmetric distribution of the wave function charge density, with its peak now again occurring in the middle 
(hfth) layer of the QDM (QDL 5 ). 

One of the most promising properties of the QDMs is the hybridization of its electron states - the 
lowest two electron states el and e 2 are bonding and anti-bonding type molecular states respectively, 
as shown in the Fig. 1 (b). This feature is a direct consequence of the strong mechanical (strain) and 
electronic couplings of the closely-spaced QD layers. Intriguingly this raises a question that as the electric 
held pushes the electron wave functions towards the edges of the QDM thereby increasing their quantum 
conhnement, does this weaken the interdot coupling, and at what held magnitude the interdot coupling 
completely breaks down? A recent experimental study [25] attempted to probe this question by measuring 
electric held dependent PL spectra. They suggested that an electrical held of 10 KV/cm maintains the 
strong interdot coupling of the QD layers, which completely break down at very large held magnitudes. 
We theoretically address this question by analysing the spatial distributions and the character of the 
electron states plotted as a function of the applied electric held magnitudes in Figs. 2 and 3. Overall 
our calculations predict that the strong interdot coupling of the QDM survives large magnitudes of the 
applied electric helds. Even at \Ep\ = 25 KV/cm and \Ea\ = 35 KV/cm, the electron states maintain 
bonding and anti-bonding type character and are still of s-type symmetry. This supports the conclusions 
from the experimental data that a strong coupling of QDLs is present in the presence of moderate electric 
held magnitudes. Our calculations also predict that the interdot coupling completely vanishes at \Ep\ = 
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FIG. 4. The interplay of the strain and the applied anti-parallel electric field is highlighted by plotting one 
dimensional (line) cuts of the lowest electron state el charge density, through the center of the QDM along the 
growth direction for the three different cases: blue line=the effect of both strain and electric field is excluded, red 
line=only the effect of strain is included and the applied electric field magnitude is still zero, green line=the effect 
of both strain and electric field is included in the simulation. In the absence of the strain and electric field, the 
bonding-like state el is symmetrically distributed along the QDM. The strain breaks the symmetry and pushes the 
state towards the bottom of the QDM. A small magnitude of the electric field {\Ea\ = 1 KV/cm) entirely cancels 
the effect of the strain and restores the symmetry of the bonding-like electron state. 

35 KV/cm and \Ea\ =45 KV/cm fields, leading to the lowest three electron states being nearly entirely 
confined in the lowest (single) QD layer of the QDM. The symmetries of the states are also analogons to 
the states of a single isolated QD layer - the lowest state el is an s-type state and the next two states (e2 
and e3) are p-type states oriented along the [-110] and [110] directions, respectively. It is worth pointing 
ont that in our calculations the difference between the Ep and Ea electric held magnitudes to achieve the 
breakdown of the interdot coupling is attributed to the asymmetric strain effect as discussed earlier in the 
section 3.1. 


3.3 — Electric field effect on the spatial confinements of the hole states 

In contrast to the electron states, which are strongly hybridized molecular states and are spread over 
multiple QD layers, the hole states are much more conhned, residing in either single QD layers or within 
the adjacent couple of QD layers, as evident from the plots of the Fig. 1 (b). The application of the electric 
helds will induce an opposing force on the hole states when compared to the electron states, due to their 
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FIG. 5. Side views of the highest four hole states (hi, h2, h3, h4) are shown in the [010]-plane as a function of the 
applied electric field parallel to the growth direction (Ep). The magnitudes of the Ep are provided in the units of 
KV/cm. The dashed lines mark the locations of the InAs wetting layers. 

opposite charge: Ep field will push the hole states towards the top and the Ea held will push the hole states 
towards the bottom of the QDM. 

Fig. 5 (a-j) plots the side views of the highest four hole states (hi, h2, h3, h4) in the [010]-plane, as a 
function of the applied electric held magnitudes in parallel to the growth direction {Ep). The dotted lines 
indicate the position of the InAs wetting layers inside the QDM geometry. Initially without the application 
of any electric held (Fig. 5 (a)), the strain ehect pushes all of the hole states towards the bottom of the 
QDM - hi is in QDL 2 , h2 is in QDL 2 and QDL 3 , h3 is in QDL 2 , and h4 is in QDL 2 and QDL 3 . Note 
that the hi and h2 states are on the right edge of the QDM, and the h3 and h4 states are on the left side 
edge of the QDM. This arrangement is related to the symmetry of the states within the QDLs and will be 
discussed later in the next section when we will describe the symmetries of the hole states with the help 
of the top view plots. This section only focuses on the conhnements of the hole states in a particular QDL 
along the growth direction. 

As we increase the magnitude of Ep held from 1 KV/cm to 35 KV/cm, the hole states move upward, 
along the growth direction of the QDM. The highest hole state hi quickly moves to the QDLy and QDLg 
layers when the held is increased from 2 KV/cm to 3 KV/cm, however it only shifts inside the topmost 
QDLg for very large helds, \Ep\ = 35 KV/cm. Our calculations also reveal that for the helds larger than 4 
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FIG. 6. Side views of the highest four hole states (hi, h2, h3, h4) are shown in the [010]-plane as a function of 
the applied electric field anti-parallel to the growth direction (Ea). The magnitudes of the Ea are provided in the 
units of KV/cm. The dashed lines mark the locations of the InAs wetting layers. 

KV/cm, all of the highest four hole states reside close to the top of the QDM. Since the ground electron 
state el for the similar field magnitudes has moved considerably towards the bottom of the QDM: at \Ep\ 
= 4 KV/cm, the el state resides in only the lowest three to four QDLs, this would considerably suppress 
the interband transition strengths involving el electron state and the hole states. We will further discuss 
this in the last section of this paper. 

Whilst the electric held Ep pushes the hole states towards the top of the QDM opposing the force applied 
by the strain helds, the held applied anti-parallel to the growth direction Ea pushes the hole states towards 
the bottom of the QDM, and therefore acts in support of the strain held. Fig. 6 (a-k) plots the side views 
of the highest four hole states (hi, h2, h3, h4) as a function of the applied Ea helds. The ground hole state 
hi, which initially resides in the QDL 2 at |F^(j|= 0 , quickly starts to shift into the lowest QDLi as the held 
magnitude increases. For 1 KV/cm < \Ea\ < 15 KV/cm, we calculate a weak hybridization of the hi state 
over the lowest two QDLs of the QDM. At \Ea\ = 25 KV/cm, the hi state is roughly equally present in 
both of the lowest QDLs (QDLi and QDL 2 ). At \Ea\ > 35 KV/cm, the hi state has completely moved 
inside the lowest QDL (QDLi). 

The shift of the hole states along the QDM growth axis not only modihes the electron-hole wave function 
overlaps and therefore interband optical transition strengths, it also changes the heavy-hole (HH) and the 


13 











light-hole (LH) character of the hole states. In the absence of electric helds, the splitting between the 
HH and LH states along the QDM growth direction is asymmetric [4]. The HH-LH splitting is minimum 
around the center of the QDM (around QDL 5 ), whereas it is much larger close to the edges of the QDM 
(inside the QDLi and QDLg). Furthermore, due to the asymmetry of the strain along the growth direction, 
the HH-LH splitting is considerably larger in the bottom QDLs compared to the top QDLs of the QDM. 
For example, in Ref. 4 the lowest layer (QDLi) exhibits an HH-LH splitting of 110 meV compare to only 
87 meV splitting in the top most QDL (QDLg). 

As we have shown that the Ep held shifts the hole states towards the bottom of the QDM, whereas 
the Ea held pushes the hole states towards the top of the QDM, this implies that the application of Ep 
held slightly increases {decreases) the HH {LH) character of the hole states, whereas the application of Ea 
held decreases {increases) the HH {LH) character of the hole states. Therefore the electric held response 
of hole states is also directional, where the HH-LH character of the hole states varies depending on the 
direction of the applied held. This change of the HH-LH character of the hole states is important for the 
polarisation-resolved optical measurements [3], as the TE mode optical strength is directly proportional 
to the HH character of the hole states, and the TM mode optical strength is directly proportional to the 
LH character of the hole states. Therefore, the application of an electric held to the QDM can engineer 
both the optical transition strengths via engineering of the electron-hole overlaps and the optical mode 
strengths via tuning of the HH-LH character of the hole states. 

3.4 — Electric field effect on the orientation of the hole states 

While the conhnement of the hole states inside the QDM directly modihes the electron-hole wave 
function overlaps and therefore the interband optical transition magnitudes, the symmetry of the electron 
and hole states is also a critical design parameter as it strongly influences the polarisation properties of 
the interband optical transitions. For a single isolated QDL, it is well-known in the literature that the 
symmetry of the lowest electron state is s-type and it is of p-type for the next (hrst) two excited states [39]. 
It was found that the strain and the atomistic conhguration of the underlying Zincblende crystals strongly 
impact the symmetry of the p-states [31]. For the lattice-mismatched InAs/GaAs QDs, the influence of the 
strong biaxial strain splits the degeneracy of the heavy-hole and the light-hole bands, lifting the heavy-hole 
band higher in energy and pushing the light-hole band lower in energy [49]. Therefore the few topmost 
hole states have dominant heavy-hole character. This is typically true for the flat QDs with low aspect 
ratios (height/base < 0.15). 

For the QDs with large aspect ratios (height/base > 0.3), the presence of a strong impact from another 
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FIG. 7. Top views of the highest four hole states (hi, h2, h3, h4) are shown in the [001]-plane as a function of 
the applied electric field parallel to the growth direction {Ep). The magnitudes of the Ep field are provided in the 
units of KV/cm. The dashed circles mark the boundaries of the InAs QD regions. 

effect was observed: the emergence of the heavy-hole pockets at the interfaces of the QDs [5, 31]. These 
heavy-hole pockets, which are formed due to a complex interaction of the strain fields and the large aspect 
ratios, result in the stronger confinement of the hole states closer to the interfaces of the QDs. Since the 
underlying strain enhances the atomistic asymmetry of the Zincblende crystals, which favours the [-110] 
direction in the plane of the QD, the heavy-hole pockets are deeper at the [-110] interfaces compared to the 
other in-plane directions. This leads to the highest few hole states being oriented along the [-110] direction. 
The large QDMs, consisting of very closely-spaced and strongly-coupled QDLs such as being studied in this 
work, could be considered as QDs with very large aspect ratios. Therefore they exhibit similar behaviour: 
the highest few hole states being conhned at the [-110] interfaces of the QDLs inside the pockets of the 
hole band edges. The plots of the top views of the highest four hole states in the [001]-plane are shown in 
Fig. 1 (c) and conhrm this notion. 

One would naively expect that the application of the electric helds in parallel or anti-parallel to the 
growth direction should not impact the orientation of the hole states in the plane of the QDLs, as it should 
primarily depend on the in-plane helds. However, as discussed in the last section, that the application of 
Ep and Ea helds results in the tunnelling of the hole wave functions from one QDL to the other QDL. 
Consequently the orientation of the hole wave functions also depend on the nature of the hole pockets 
belonging to that particular QDL in which a hole state resides at an applied held magnitude. Therefore 
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FIG. 8. Top view of the highest four hole states (hi, h2, h3, h4) are shown in the [001]-plane as a function of 
the applied electric field anti-parallel to the growth direction (Ea). The magnitudes of the Ea are provided in the 
units of KV/cm. The dashed circles mark the boundaries of the InAs QD regions. 

changes in the orientation of the hole states could be observed as a function of the electric helds. 

To investigate the effect of electric helds on the orientations (symmetries) of the hole states, we vary Ep 
and Ea helds and plot the top views of the held-dependent hole wave function in the [OOlj-plane as shown 
in Figs. 7 (a-j) and 8 (a-k). The dotted circles roughly mark the boundaries of the QDLs. The topmost 
row in Fig. 7 shows that the application of an electric held in parallel to the growth direction, Ep, does 
not impact the orientation of the ground hole state hi. Therefore the hi state remains oriented along the 
[-110] direction for all of the applied held magnitudes. Interestingly the other three hole states h2, h3, 
and h4 exhibit some changes in their orientations when the magnitude of the Ep held is increased. For 
example, the h2 state is oriented along the [-110] direction for \Ep\ < 3 KV/cm and is oriented along the 
[110] direction for \Ep\ > 3 KV/cm. Even a stronger variation in the orientations of the h3 and h4 states is 
observed. At 35 KV/cm, the interdot coupling breaks down and the electron and hole states are conhned 
in the individual QDLs. Therefore, the orientation of the hole states for 35 KV/cm is what we expect for 
a single QD: the hole states are alternatively oriented along the [-110] and the [110] directions [50]. This 
implies that the in-plane polarisation modes TEpioj and TE[_iio] will be of equal magnitudes, despite the 
conhnements of the hole states being strongly directional, as previously reported for the single QDs [50]. 

For the highest hole state hi, the behaviour is same for the held applied anti-parallel to the growth 
direction, Ea'. the hi state is oriented along the [- 110 ] direction independent of the applied held magnitudes. 
Therefore we conclude that for the hi state, the orientation is independent of both, the held magnitude and 
the held direction. Consequently if the optical transition is designed to involve only the hi state, the only 
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dominant optical mode would be the TE[_iio] mode. The h2 state exhibits a slightly different behaviour 
under the influence of Ea held, when compared to the Ep held: the h2 state is also oriented along the [-110] 
direction for all of the applied Ea magnitudes. However, it switches conhnement in the HH pockets from 
upper-left corner of the QDL to the lower-right corner for the certain values of the helds as shown in the 
second row of the Fig. 8. Finally we also calculate that the h3 and h4 states also exhibit slightly diherent 
Ea held response, compared to the Ep held. One important diherence is the emergence of p-type hole h3 
and h4 states at \Ea\ =45 KV/cm held. 

Overall, from the examination of the held dependent hole state orientations, we conclude this section 
by stating that the helds applied along the perpendicular directions to the growth direction modify the 
symmetries of the hole states in the plane of the QDM. This unintuitive result is attributed to the complex 
behaviour of the HH pockets found at the interfaces of the QDLs inside the QDM. 

3.5 — Transition Energies and Band-gap wavelengths 

The electric held dependence of the ground state transition energy (Transe=el-hl) is plotted in Fig. 9 
(a), where the positive and negative values of the electric held indicate the Ep and Ea helds, respectively. 
For a QD system, the Trans^ value is represented by the following equation [29]: 

TranSe{\Ep/a\) = TranSe{0) + p\Ep/a\ -F /3\Ep/a\‘^ (3) 

where TranSe{\Ep/a\) is the transition energy value at \Ep/a\ held, TranSe{0) is the transition energy 
at reference held (here we assume it to be zero), p depends on the built in dipole moment due to the 
electron-hole separation induced by the applied electric helds, and /? is a measure of the polarizibility of 
the electron and hole wave functions. For a typical single InAs QD, the Transe dependence on the electric 
held is parabolic [29, 51], indicating a large magnitude for the polarizibility f3 and a very small magnitude 
for the dipole moment p - a direct consequence of the small space available for the electron and hole 
separation. For QD stacks, made up of a large number of QD layers, electric helds are capable of pushing 
electron and hole wave functions far away from each other (even at the opposite ends of the stack for large 
helds), and therefore can lead to large separations between the electron and hole wave functions. This 
results in much larger values for the dipole moment p compared to the polarizilibilites (3 - hence the held 
response of the QD stacks is characterized by a linear dependence of the TranSe on the applied electric 
helds. 

Fig. 9 (a) overall demonstrates a linear dependence of the Transe value on the applied held magnitudes 
for the studied QDM. For the large values of the applied electric helds, \Ep/a\ > 2 KV/cm, the electric 
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field dependence of the Transe is nearly perfectly linear and therefore can be represented by the following 
equation: 


TranSe{\Ep/a\) = TranSe{0) + p\Ep/a\ (4) 

The black dotted lines in the Fig. 9 (a) plots the hts from the equation 4, with the values of p used as -3.65 
meV.cm/KV and -3.75 meV.cm/KV for the positive and the negative electric helds, respectively. These 
values of the p are signihcantly larger compared to the single QD values - several orders of magnitude 
smaller values have been reported in the literature for the single QDs [29]. 

It is also noted that for the small values of the applied electric helds, |i?p/a| < 2 KV/cm, the TranSe 
values exhibit a parabolic dependence on the electric held magnitudes. This is clear from the inset plot 
in Fig. 9 (a), where the dotted black line shows the htting from the equation 3 with the values of p = 
-3.75 peV.cm/KV and (3 = -0.25 meV.cm^/KV^. While the magnitude of p is decreased by three orders of 
magnitude, the value of the polarizibility (3 is signihcantly large. This is also considerably larger compared 
to the measured values of -2 to -5 peV.cm^/KV^ for the single QDs [51]. Therefore we conclude that a 
strongly coupled QDM can be operated in both regimes: a large polarizibility {dipole moment) regime for 
the small {large) magnitudes of the applied electric helds. 

It should also be noted from the inset of the Fig. 9 (a) that the maximum value of the Trans^ parameter 
occurs at 2 KV/cm electric held, rather than 0. This slight asymmetry is a direct consequence of the 
underlying strain held that pushes the electron state el towards the top of the QDM and the hole state hi 
towards the bottom of the QDM as discussed earlier in the section 3.1 (see also Fig. 1 for reference). 

Another parameter of interest for the optoelectronic devices is the band-gap wavelength. Typically 
large strain helds blue shift the optical wavelengths of the QD systems and therefore make them unsuitable 
to operate at the optical telecommunication wavelengths (1.3 pm and 1.55 pm). It has been a great 
challenge to achieve the optical emissions from the QD systems at 1.3 pm and 1.55 pm wavelengths, 
where several diherent techniques have been used in the literature to overcome this problem. The most 
widely used QD system ohering the long wavelength emission is QD inside a quantum well structure, where 
strain engineering results in the red shift of the PL peak [49, 50]. The application of the electric helds also 
induces large shifts in the electron and hole conhnement energies, which is commonly known as the quantum 
conhned Stark effect (QCSE). The QCSE results in reducing the band-gap energy of the QD system and 
therefore red shifts the band-gap wavelengths. Here we investigate the band-gap wavelength increase as 
a function of the applied electric helds. For a single QD, the QCSE is only dependent on the bending of 
the band edges and therefore it primarily exhibits a linear dependence on the magnitude of the electric 
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FIG. 9. (a) Transition energies (el-hl) are plotted as a function of the magnitudes of the applied electric fields. For 
the large electric field magnitudes iFlp/al > 2 KV/cm, the held dependence of the transition energy is nearly linear 
indicating negligible polarizibility (/5). For < 2 KV/cm, the transition energy exhibits a perfect parabolic 

dependence on the electric held implying very small dipole moment (p) (see the inset), (b, c) Band-gap wavelength 
(1240/el-hl) in nm is plotted as a function of the magnitudes of the applied electric helds: (b) the held direction 
is parallel to the growth direction (Ep) and (c) the held direction is anti-parallel to the growth direction {Ea). The 
inset plots highlight the wavelength dependence for the small magnitudes of the electric helds. 


field [28]. Similar behaviour has been observed for weakly-coupled bilayer QD systems [30]. However for 
a strongly-coupled QDM, such as the one being investigated in this study, the strong interdot couplings 
lead to a non-trivial dependence of the QCSE on the magnitude of the electric helds. We have previously 
discussed in the section 3.1 that the electron states are strongly coupled hybridized states whereas the hole 
states are much more conhned and typically present only in one or two QDLs of the QDM. This implies 
that the electric held would lead to the tunnelling of the hole states between the QDLs and therefore adds 
an extra level of complexity in the QCSE behaviour. 

We plot the band-gap wavelengths for the QDM as a function of both Ep and Ea electric held magnitudes 
in Figs. 9 (b) and (c), respectively. The band-gap wavelengths are calculated by converting the energy 
diherence of the lowest electron state el and the highest hole state hi to a corresponding wavelength on 
the nm scale: (1240/el-hl). Since the strain distribution along the growth direction is asymmetric, the 
electric held induced QCSE in the band-gap wavelengths is also to be diherent for the Ep and Ea helds. 
The diherence mainly arises from the tunnelling of the hole states as a function of the electric helds. For 
the Ep held, the hole state hi shifts from the bottom to the top of the QDM and therefore experiences a 
strong interdot tunnelling ehect. This is rehected in the inset of Fig. 9 (a), where a strong deviation from 
the typical linear increase of the band-gap wavelength is calculated. However this ehect only lasts for very 
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small magnitudes of the Ep field. At 2 KV/cm field, the hole state hi has already moved to the top of 
the QDM, so the tunnelling related effect disappears. Further increase in the Ep field magnitude leads to 
roughly linear increase in the band-gap energy. 

The application of the field anti-parallel to the growth direction, Ea, does not show much deviation from 
the linear behaviour, even for the small magnitudes of the field. This is because the hi state is already 
close to the base of the QDM at zero field, and therefore it only experiences a weak tunnelling effect when 
\Ea\ is switched on. The inset of the Fig. 9 (b) demonstrates this behaviour for \Ea\ < 2 KV/cm. 

3.6 — Interband optical transition strengths 

For a variety of optoelectronic devices such as semiconductor optical amplifiers (SOAs), polarisation- 
resolved interband optical transitions are a critical design component. Several experimental [3, 24] and 
theoretical [4, 52] studies have explored this parameter to attain isotropic polarisation response, where the 
optical mode strengths are independent of the polarisation direction. Polarisation-dependent interband 
optical transition strengths (TE and TM modes) depend on two factors: (1) the spatial overlap between 
the electron and the hole wave functions and (2) the heavy-hole (HH) and the light-hole (LH) character 
of the hole states. In QDMs, the strong coupling between the QDLs lead to enhanced LH character of 
the valence band states. Furthermore the spatial overlap between the electron and hole states is decreased 
due to the interfacial confinements of the hole states (see Fig. ?? for reference). The two factors lead 
to strong modifications of the TE and TM polarisation modes. In particular, the previous experimental 
measurements [3] and theoretical calculations [4] have revealed a unique property for the QDMs, where the 
TM mode strength becomes comparable to the TE mode strength - a characteristic critically important 
for the design of semiconductor optical amplihers (SOAs). This property is not accessible from single QDs, 
where the TE mode is always dominant and the TM mode is very weak, resulting in a strong anisotropy 
of the polarized absorption spectra [50]. 

In this section, we examine the polarisation resolved interband transition modes as a function of the 
applied electric helds. Since the electric held has been shown to strongly modify the electron and hole spatial 
conhnements, the optical mode strengths are also expected to exhibit correspondingly large variations. For 
a complete spectra of the optical modes, we investigate the TE modes in the four (high symmetry) in¬ 
plane directions, which are of particular relevance to the experimentalist: [-110], [110], [100], and [010]. The 
corresponding TE modes are labelled with the polarisation directions in the subscripts. The TM mode 
is always calculated along the growth [001] direction of the QDM. Fig. 10 plots the polarisation mode 
strengths as a function of the applied electric fields: (a) when the field is parallel to the growth direction. 
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FIG. 10. Interband optical transition modes (TE and TM modes) calculated form the equations 1 and 2 are plotted 
as a function of the applied electric fields: (a) The applied field direction is parallel to the growth direction of the 
QDM and (b) the applied field direction is anti-parallel to the growth direction of the QDM. 


\Ep\ and (b) when the held is anti-parallel to the growth direction, \Ea\- 

The application of the electric helds shifts the electron and hole wave fnnctions in the opposite directions 
and therefore reduces the electron-hole wave function overlaps. Therefore, overall all of the optical mode 
strengths decrease as a function of the electric held. A careful analysis of the plots of Fig. 10 highlights the 
following insights: (1) The decease in the optical mode strengths for small magnitudes of the electric helds 
is smoother for the Ea held, when compared to the Ep held. This is simply due to the fact that the hole 
wave functions experience signihcantly weaker interdot tunnelling under the application of the Ea held, 
compared to the Ep held. (2) The TE[_iio], TEjioo], and TEpio] optical modes remain much stronger than 
the TM[ooi] mode for all magnitudes of the applied helds, and it is also independent of the held direction. 
Therefore, isotropic polarisation is not accessible in these polarisation directions with or without helds. 

(3) The magnitudes of the TE[ioo] and TE[oio] modes are always same, independent of the electric held 
magnitudes and directions. This is due to the fact that the electric helds in this study are only applied 
along the growth [001] direction and therefore do not have any directional ehect for the in-plane polarisation 
modes. Furthermore, the hole wave functions are conhned along the [110] or [-110] directions, therefore the 
[100] and the [010] directions are equivalent with respect to the electron and hole wave function overlaps. 

(4) At the zero held value, TEjno] > TM[ooi]. However for small magnitudes of the electric helds (for both 
Ep and Ea), approximately varying from 2 KV/cm to 4 KV/cm, the TEpioj mode strength is roughly equal 
to the TM[ooi] mode strength. This implies that the application of the small electric helds can be used to 
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realize isotropic polarisation from QDMs. (5) Due to the large separation between the electron and hole 
wave functions, the optical mode strengths rapidly decrease as a function of the applied electric helds. At 
15 KV/cm, we calculate that the optical mode strengths have been reduced by approximately two to three 
orders of magnitudes. Therefore, very weak interband optical absorption is expected at the large electric 
helds. This is expected to be important for the design of photovoltaic devices such as intermediate-band 
solar cells where optically generated electron-hole pairs will be prohibited from reabsorption via interband 
transitions. 

4 - CONCLUSIONS 

In conclusion, we have provided a comprehensive understanding of the electric held control of the elec¬ 
tronic and optical properties of the multi-layer strongly-coupled quantum dot molecules (QDM). Electric 
helds are applied in both parallel {Ep) and anti-parallel {Ea) directions with respect to the growth axis 
([001]) of the QDM. It is found that the asymmetric ehect of the strain on the electronic wave functions can 
be entirely balanced by the application of a small electric held of the order of 1 KV/cm along a direction 
anti-parallel to the growth direction. We have provided a detailed investigation of the electron and hole 
wave functions under the inhuence of electric helds by plotting the top and side views, which would serve 
as a useful guide for the experimentalists working on the devices based on the QDMs. Our calculations 
highlight a complex interplay between the applied electric helds and the strong interdot couplings, which 
control the spatial conhnements of the electron and hole wave functions in diherent QDLs of the QDM. It 
is revealed that the interdot couplings survive large electric held magnitudes and a complete break down 
occurs at \Ep\ >35 KV/cm or |Ea| < 45 KV/cm. At these large electric helds, the electronic states reduce 
to single quantum dot states conhned at the opposite ends of the QDM. This mimics a transformation from 
a type-I band structure to a type-II band structure, which is a requirement for the design of intermediate- 
band solar cell devices. By analysing the ground state transition energy data, we highlight an interesting 
property of the QDMs that it could be operated in two regimes: for the small electric held magnitudes, the 
QDM exhibits large polarizibilities and negligible dipole moment strengths; in contrast for the large electric 
held magnitudes, the QDM demonstrates negligible polarizibilities and the held response is dominant by 
the dipole moment. The band-gap wavelength of the QDM red shifts as a function of the applied electric 
helds due to the quantum conhned Stark ehect (QCSE) and increases to 1.3 fj,m at the 15 KV/cm held. 
However, the reduced overlap of the electron-hole wave functions results in the optical transition strengths 
reduced by two to three orders of magnitude, thereby a much smaller optical gain/oscillator strength is 
expected under the inhuence of the electric helds. We have also provided a detailed investigation of the 
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polarisation-resolved interband optical transition modes (TE and TM modes) as a function of the applied 
electric fields. Overall the polarisation mode strengths decrease as a function of the applied electric helds 
due to the reduced electron-hole wave function overlaps. However more importantly, the application of the 
small electric helds is found to be helpful to realize the isotropic polarisation response from the QDMs, 
with TEjiio] ~ TM[ooi]. This condition is crucial for the design of semiconductor optical amplihers (SOAs) 
where the isotropic polarization is highly desirable. Overall, our multi-million atom calculations provide a 
critical and reliable theoretical guidance for the experimentalists on a topic where little theoretical knowl¬ 
edge is available in the existing literature. The physical insights presented will be highly useful to engineer 
electronic and optical properties of the QDMs for the design of a variety of optoelectronic, photovoltaic, 
and quantum information science devices. 
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